Net baryon density in Au + Au collisions at the Relativistic Heavy Ion Collider 
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We calculate the net baryon rapidity distribution in Au+Au collisions at the Relativistic Heavy 
Ion Collider (RHIC) in the framework of the Parton Cascade Model (PCM). Parton rescattering 
and fragmentation leads to a substantial increase in the net baryon density at mid-rapidity over the 
density produced by initial primary parton-parton scatterings. The PCM is able to describe the 
measured net baryon density at RHIC. 



First experiments at the Relativistic Heavy Ion Col- 
lider (RHIC) have shown that the matter created in 
the central rapidity region contains a significant excess 
of baryons over antibaryons. While a slight baryon 
excess was not unexpected, the magnitude of the net 
baryon multiplicity density dN^_-g/du « 19 ± 2 P, 
at Vi^jv = 130 GeV and « 14 ± 4 jg at y/sj^j^ = 200 
GeV is higher than what many theoretical models had 
predicted jj]. In particular models in which the deposi- 
tion of energy at midrapidity is driven by quasiclassical 
glue fields [J or fragmenting color flux tubes 0], which 
produce quarks and anti-quarks in equal abundance, un- 
derpredicted the data. On the other hand, models invok- 
ing baryon junctions 0, HI for the transport of baryon 
number from the beam rapidity into the central region 
y have done remarkably well. The baryon junc- 
tion mechanism was originally proposed as a means to 
understand baryon number annihilation and stopping in 
elementary p + p and p + p reactions 9] . In this let- 
ter we shall address the question whether only baryon 
junctions provide a mechanism capable of explaining the 
RHIC data, or whether they can also be understood in 
the framework of a more conventional picture, based on 
parton distributions and pQCD driven multiple interac- 
tions. 

Two effects based on established physics can contribute 
to the baryon excess at midrapidity. First, the mea- 
sured parton distribution functions in the nucleon are 
well known to exhibit a substantial asymmetry between 
quark and antiquark distributions at moderately small 
values of Bjorken-x {x w 0.01) relevant for RHIC energies 
( Fig.^ see below). Second, it is known from experiments 
with p -\- A collisions that multiple scattering is quite ef- 
fective in transporting baryons to smaller rapidities |lfll | . 
The parton cascade model (PCM) [HH 0| provides 
a natural framework for exploring the consequences of 
these two effects quantitatively. 

Devised as a description of the early, pre-equilibrium 
phase of a nucleus- nucleus collision at relativistic energy, 
the PCM does not include a description of hadronization 
and the subsequent scattering among hadrons. These 
late-stage processes, however, are not expected to alter 



the distribution of net baryon number with respect to 
rapidity, because baryon diffusion in a hadronic gas is 
slow |l4| and net baryon number is locally conserved. 
We therefore believe that these limitations of the PCM 
approach should not stand in the way of an adequate ex- 
planation of the net baryon distribution. In this work, we 
present calculations of the net baryon multiplicity den- 
sity distribution in Au+Au collisions at 130 and 200 GeV 
center-of-mass energy per nucleon pair. 

The fundamental assumption underlying the PCM is 
that the state of the dense partonic system can be 
characterized by a set of one-body distribution func- 
tions Fi{x^,p°'), where i denotes the flavor index [i ~ 
g,u,u,d,d, . . .) and x^,p°' are coordinates in the eight- 
dimensional phase space. The partons are assumed to 
be on their mass shell, except before the flrst scatter- 
ing. In our numerical implementation, the GRV-HO 
parametrization [isj is used, and the parton distribu- 
tion functions are sampled at an initialization scale 
to create a discrete set of particles. Partons generally 
propagate on-shell and along straight-line trajectories be- 
tween interactions. Before their first collision, partons 
may have a space-like four-momentum, especially if they 
are assigned an "intrinsic" transverse momentum. 

The time-evolution of the parton distribution is gov- 
erned by a relativistic Boltzmann equation: 

P'-^,F,{^^P)=CAF] (1) 

where the collision term Ci is a nonlinear functional of 
the phase-space distribution function. The calculations 
discussed below include all lowest-order QCD scattering 
processes between massless quarks and gluons . A low 
momentum-transfer cut-off p™'" is needed to regularize 
the infrared divergence of the perturbative parton-parton 
cross sections. Additionally, we include the branchings 
q ^ qg, q ^ qj, g ^ gg and g ^ qq [i3|- The soft 
and coUinear singularities in the showers are avoided by 
terminating the branchings when the virtuality of the 
time-like partons drops below fiQ — 1 GeV. Some of these 
aspects were originally discussed in [l^. The present 
work is based on our thoroughly revised, corrected, and 
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FIG. 1: Net baryon content of the partonic distribution func- 
tion of gold-nucleus at the factorization scales (Qo) of 0.50 
GeV^ (sohd curve), 1.00 GeV^ (dashed curve) and 2.00 GeV^ 
(dot-dashed curve), for the GRV-HO Il5^ parametrization. 
Parton shadowing is not included. 



extensively tested implementation of the parton cascade 
model, named VNI/BMS 0. 

Figure ^ shows the net baryon-number distribution 
X A[fq{x) — f-q{x)] /3 with the Bjorken-x variable. We 
have defined, 



(2) 



where Fi{x) is the iso-spin averaged parton-distribution 
function for a nucleon in the gold nucleus, so that 



P^{^) = ^Ff{x) + 



{A_Z) 
A 



(3) 



where the superscripts P and N stand for protons and 
neutrons. Thus the figure brings out the asymmetry be- 
tween the quark and antiquark distributions at moder- 
ately small values of Bjorken-x {x k, 0.01) relevant for 
energy deposition at midrapidity at RHIC energies. The 
various lines denote different initialization scales for the 
parton distribution functions. This provides the basis 
for the evaluation of the dynamical evolution of the stop- 
ping through parton collisions, sometimes accompanied 
by gluon radiation and quark pair creation. 

The partons in a fast moving nucleon or a nucleus are 
distributed both in the transverse and longitudinal direc- 
tions. Due to the large longitudinal Lorentz contraction, 
the longitudinal momentum pz = xP is the most interest- 
ing variable. Before or immediately after the first colli- 
sion of a parton, one can relate Pz to the rapidity variable 
y = Y + lnx + \n{M/Qs), where Y is the rapidity of the 
fast moving nucleon, M is the nucleon mass, and Qs de- 
notes the typical transverse momentum scale. Depending 
on the picture of the initial state, Qs is either given by 



the average intrinsic virtuality, often called the satura- 
tion scale |l8j. or by the typical transverse momentum 
given to the parton in the first interaction which brings 
it onto the mass shell. In any case, |ln(A//(5s)| < 1/2 
for Au-|-Au collisions at RHIC. 

If the partons in the two gold nuclei were to decohere 
completely upon passing through each other without any 
further interaction, the resulting net baryon rapidity dis- 
tribution would be approximately given by figure^- pre- 
dicting a contribution to the net baryon density of about 
7 at mid-rapidity (with Y = 5.4 for y/sj^j^ =200 GeV) 
due to each nucleus. This is seen as follows. The baryon 
number in the nucleus A is given by 



[X fq{x) - X fqix)] dlliX 



(4) 



so that the net baryon-distribution "contained" in the 
nucleus is 



dN 



1 



dy J B-B 3 



XA [fq{x)-fq{x)] 



(5) 



We shall see that the parton-interactions incorporated in 
the PCM remain consistent with this intuitive picture. 

A crucial parameter of the PCM is the low momentum 
transfer cut-off p™'". Under certain assumptions this pa- 
rameter can be determined from experimental data for 
elementary hadron-hadron collisions [l^ Eol l2ll | . In the 
environment of a heavy-ion collision, color screening will 
destroy the association of partons to particular hadrons, 
since the density of free color charges is so high that the 
color screening length becomes smaller than the typical 
hadronic scale. In a previous publication we have 
established a consistency limit for the allowed range of 
p™'" by calculating the screening mass hd the pro- 
duced parton matter as a function of the cut-off p™'" and 
demanding (p™'" > ^d)- The initialization scale and low 
momentum cut-off of the pQCD cross sections were cho- 
sen as Ql = (p??'")2 = 0.50 GeV2 for ^i^j^ = 130 GeV 
and 0.59 GeV^ for ^i^j^ = 200 GeV, respectively The 
energy scaling of this parameter agrees with that deter- 
mined by Eskola et al. for the geometric minijet 
saturation model. 

Figure|21shows the PCM predictions for the net baryon 
rapidity distributions for Vs^jv ~ -'^30 GeV (upper 
panel) and 200 GeV (lower panel), respectively. The ini- 
tial projectile and target rapidities are ±4.9 for Vs^vat = 
130 GeV and ±5.4 for Vsatat = 200 GeV. The calcu- 
lations are done without assigning any intrinsic kx to 
the partons - initially all partons move with the velocity 
[3 — ±P^ I Ea of the nucleus in the center-of-mass frame, 
where Pz is the momentum and Ea is the energy of the 
nucleus. 

Crosses in Fig. El denote a calculation in which the 
PCM has been restricted to primary-primary parton 
scatterings, and therefore reflects a calculation in which 
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FIG. 2: 

tions at 



net baryon rapidity distributions for Au+Au reac- 



130 GeV (top) and 



200 GeV bot- 



tom. Crosses denote a calculation in which the PCM has 
been restricted to primary parton scatterings, diamonds in- 
clude parton rescattering and squares include rescattering and 
parton fragmentation. The solid lines show the net baryon 
content of the partonic distribution functions for gold nuclei, 
scaled by an average liberation factor of 0.4 . The band at 
ycAi shows the range of experimental estimates for the net- 
baryon density by STAR, BRAHMS and PHENIX. 



each parton is allowed to scatter only once. Already 
one hard collision is sufficient to deposit a net surplus 
of quarks into the mid-rapidity region, resulting in a net 
baryon density at ycm = of 6.3 for 130 GeV and 5.0 at 
200 GeV. For comparison, the net baryon number distri- 
bution for each colliding nucleus, scaled by a factor 0.4 
from the distribution shown in Fig.^ is shown as a solid 
line. The remarkable agreement demonstrates that the 
net baryon number distribution produced by first parton- 
parton collisions is predetermined by the initial parton 
structure of the nuclei. The factor 0.4 is the average "hb- 
eration factor" c for partons in the PCM for the selected 
parameters. This factor is consistent with predictions by 
some gluon saturation models |2^ . 



The diamonds in Fig.[21represent a calculation with full 
parton-parton rescattering. Allowing for multiple parton 
collisions increases the net baryon density at mid-rapidity 
roughly by a factor of two at 130 GeV and by 75% at 
200 GeV, filling up the dip around mid-rapidity. This 
trend continues when parton fragmentation is included 
(squares): at 130 GeV fragmentation processes add an- 
other 50% to the net baryon density at mid-rapidity, 
bringing it up to about 18, whereas at 200 GeV the net 
baryon density increases to near 14. The rapidity change 
of a quark in each subsequent collision after its liberation 
in the first hard scattering yields an average rapidity shift 
of roughly 0.65 units per collision. 

The band around mid-rapidity in Fig. |21 denotes the 
range of experimental estimates for the net-baryon den- 
sity at mid-rapidity for 130 GeV and 200 GeV, respec- 
tively P, 0, 13 ■ These estimates depend on how the ex- 
trapolations from {p—p) and (A— A) to (B—B) have been 
carried out, the extreme scenarios being the assump- 
tion of full isospin equilibration for baryons, {B — B) = 
2 X (p—p), vs. the full inclusion of the isospin asymmetry 
of the initial state, [B — B) = A/Z x {p — p)- Note that 
our calculation does not suffer from this uncertainty. 

Figure O shows the net baryon density at mid-rapidity 
as a function of the momentum cut-off p™'" for Au-|-Au 
collisions at -s/stvw — ^00 GeV. Again, crosses denote the 
calculation restricted to primary-primary collisions, dia- 
monds represent the full rescattering mode and squares 
include the effect of parton fragmentations. The observed 
power law dependence of the net baryon density as a 
function of Qo stems from the properties of the pQCD 
cross sections in the PCM. The absence of a saturation at 
small values of Qq indicates that not all valence quarks 
are "liberated" in the range of cut-off values considered 
here. Indeed, we find that the liberation factor for quarks 
in the nuclear parton distributions varies from about 0.7 
for X > 0.1 to about 0.2 for x « 0.01. Figure can be 
used to rescale the PCM prediction for the net baryon ra- 
pidity distribution in Fig. |21 (bottom part) to other values 
of Qo- For the range of Qo values extracted from RHIC 
data the PCM is well able to describe the measured net 
baryon excess. 

This success raises the obvious question: what does our 
model predict for the net baryon distribution at lower en- 
ergies? We have reported |2J| recently that the partonic 
cascades provide only a dilute medium at SPS energies 
iV^NN ~ 20 GeV), which does not support enough mul- 
tiple scattering among partons to justify a perturbative 
treatment. This implies that partonic cascades having 
Pt > Pt™ constitute only a small part of the dynamics 
of the collision at SPS energies. Indeed we find that the 
net-baryon multiplicity at central rapidity in Pb-|-Pb col- 
lision at -s/sTVAf = GeV is only about 20% of the ex- 
perimental value estimated by the NA49 experiment [i^ . 
Similar considerations also apply to p+p collisions, which 
do not produce a dense partonic medium where color 
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FIG. 3: Initialization scale and cut-off dependence of the 
net baryon density at mid-rapidity for Au+Au collisions at 
yi^^ = 200 GeV. 

screening occurs at a perturbative scale p^^"" > Aqcd- 
The transport of net baryon number then must be due to 
nonperturbative mechanisms. The situation is dramat- 
ically different in Au + Au collisions at RHIC energies, 
where multiple parton scattering at px > p™'" produces a 
medium, in which color is screened at a sufficiently short 
distance to allow for a choice of p™'" in the domain of 
pQCD. 

In conclusion, we find that the parton cascade model 
predicts a net baryon excess at mid-rapidity in Au-|-Au 
collisions at RHIC, which is in qualitative agreement with 
the measured values. Two mechanisms are driving this 
excess: One is the presence of a net baryon density in 
the initial state parton distributions at Bjorken-x around 
0.01 reflecting the size of the valence quark component in 
this range of x. The other important factor is the rescat- 
tering among partons, which transports more partons, 
and hence additional net baryon number, to mid-rapidity. 
This transport mechanism increases the net baryon num- 
ber density well into the range of measured values. 
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